American Fern Journal 92(3): 185—213 (2002) 


Differentiation of Eastern North American 
Athyrium filix-femina Taxa: Evidence 

From Allozymes and Spores 

Carol Lynn Kelloff 1 and Judith Skog 

Department of Biology. George Mason University, Fairfax, VA 22030 
US National Herbarium. Smithsonian Institution, NHB 166. Washington, DC 20560 

Laura Adamkewicz 

Department of Biology, George Mason University, Fairfax, VA 22030 

Charles R. Werth 

Department of Biological Sciences. Texas Tech University, Lubbock, TX 79409 


Abstract. — Athyrium filix-femina (Lady Fern) comprises a complex of homoploid (n = 40) taxa, 
distributed over much of the northern hemisphere and extending into South America, whose 
evolutionary relationships are poorly understood and whose taxonomic treatment is problematic. 
The A. filix-femina complex of North America comprises as many as four taxa with overlapping 
ranges and provides an especially suitable context for exploring patterns and processes of diver¬ 
gent evolution and its taxonomic consequences in ferns. We addressed differentiation of two 
eastern North American taxa distinguished on the basis of growth form, frond shape, and spore 
color, and most recently treated as varieties A. angustum and A. asplenioides (Northern and 
Southern Lady Fern respectively). Although, these two taxa have been long perceived as closelv 
related, they have been known to intergrade and recombine to form a hybrid zone in their rela¬ 
tively narrow region of overlap. This perception is supported by the data from the present study. 
Gollections from 17 populations, 9 of A. angustum from Quebec to Pennsylvania and 8 of 
A. asplenioides from New Jersey to North Garolina. were examined for spores using LM, 
SEM, and TEM, and/or allozymes (16 loci coding 10 enzymes). The two taxa exhibited highly 
distinct perispore surfaces: all A. angustum individuals had papillose surfaces, whereas most 
A. asplenioides individuals were rugose with a reticulum of inflated folds. Spores from the 
northernmost A. asplenioides population sampled (Shirley, NJ) showed varying degrees of inter¬ 
mediacy suggestive of introgressive hybridization with A. angustum. Levels of allozyme poly¬ 
morphism in populations (means: P=36.5%, A=1.97, H t = 0.129) were near means for 
angiosperms and ferns. Genotype frequencies at most loci in all populations were in Hardv- 
Weinberg equilibrium indicating an outcrossing mating system. Most alleles were shared among 
all populations. However, at the four most polymorphic loci ( Idh-1. Pgi-2. Pgm-2. and Tpi-2 ) 
allele frequencies were significantly divergent between populations of A. angustum and A. asple- 
nioides, especially ldh-l which approached fixation for alternate alleles. Values for F ST ranged 
from 0.008 to 0.459 for individual loci (0.255 across loci) with especially high values for Idh-1, 
Pgi-2, Pgm-2 , and Tpi-2. Hierarchical F ST analysis indicated that differences between the two 
taxa (F X y = 0.216) accounted for most allele frequency divergence among populations 
(F xy = 0.2 38). UPGMA analysis of paired Rogers’ Similarity (S) values resulted in two principal 
clusters each comprising populations of one taxon. Populations of A. angustum and A. asple¬ 
nioides were joined within their clusters at S=0.938 and S=0.945 respectively, while the two 
taxon clusters were joined at S = 0.848. The spore and isozyme data indicate substantial diver¬ 
gence between A. angustum and A. asplenioides. suggesting that they merit distinction at the 
rank of subspecies or species. Additional study of populations in their region of sympatrv is re¬ 
quired to determine the nature and extent of hybridization. 
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Athvrium filix-feminci (L.) Roth (Lady Fern) comprises a wide-ranging 
complex of divergent, homoploid (n = 40) taxa of allopatric or parapatric 
distribution, variously divided and treated at different ranks (species, sub¬ 
species, or variety) bv different authors. At least three taxa occur in North 
America, which Butters (1917) treated as separate species "amply distinct 
from each other,” although he considered A. filix-feminci from western North 
America as conspecific with the "true" A. filix-feminci of Europe. Wherry 
(1961) followed Butters (1917) in recognizing as distinct species the western 
A. filix-feminci (specifying var. sitchense Ruprecht), the northeastern 
A. angustum (Willd.) Presl., and the southeastern A. asplenioides (Michx.) 
A. Eaton, although noting that they "intergrade to such an extent as to defy 
any simple classification." Lellinger (1985) treated these taxa as subspecies 
of a single globally distributed species, i.e., western A. filix-feminci ssp. cvclo- 
sorum (Rupr.) C. Chr., northeastern A. filix-feminci ssp. angustum (Willd.) 
Clausen and the southeastern A. filix-femina ssp. asplenioides (Michx.) 
Hulten. More recently, Kato (1993), stating that "the delimitation and infra- 
specific classification of A. filix-femina need detailed study," treated these 
taxa at varietal rank, recognizing four North American varieties: in the west, 
northern A. filix-femina var. cyclosorum Ruprecht (= var. sitchense), a more 
southerly A. filix-femina var. californicum Butters (which Butters, 1917 had 
noted as distinctive but treated as a variety of A. filix-femina), and in the 
east, northern A. filix-femina var. angustum (Willd.) G. Lawson, as well as a 
southern A. filix-femina var. asplenioides (Michx.) Farwell. 

These differences in classification reflect neither strongly divergent view¬ 
points among authors as to the nature of species versus infraspecific taxa, nor 
differing insights following acquisition of compelling data. Rather, they reflect 
the best judgment of individual authors in the face of a lack of critical informa¬ 
tion as to character consistency, degree of intergradation, and genetic relation¬ 
ships among taxa in a group noted for a high degree of variability (Schneller 
and Schmid, 1982). As a contribution toward elucidating the nature of the taxa 
composing the A. filix-femina complex, the present investigation addressed 
the distinctness of the two parapatric taxa of eastern North America, A. filix- 
femina var. angustum and A. filix-femina var. asplenioides (referred to hence¬ 
forth informally as bare epithets angustum and asplenioides respectively). 
Characters differentiating these taxa. provided by Butters (1917) and by and 
large echoed in more recent treatments (Lellinger, 1985; Kato. 1993). include 
rhizome habit, leaf shape, and notably color and surface features of the spores 
(Table 1). The goal of the present study was to evaluate more fully the degree 
to which these putative taxa are distinct by 1) characterizing fine features of the 
spores using electron microscopy, and 2) surveying allozyme variation across 
a north-to-south transect through the eastern portion of the range of both taxa. 


Materials and Methods 

Collections. —All observations were made on plants collected by CLK and 
CRW from 17 localities, nine A. angustum and eight A. asplenioides (Table 2). 



Table 1 . Comparison ol morphological Features of the two putative eastern North American taxa in the Athyrium filix-femina complex. Compiled 
from taxon descriptions in Butters (1917), Lellinger (1985), and Kato (1993), 


Character 

Rhizome orientation 
Scales: 


Stipe length relative to 
lamina 

Lamina shape 
Frond base 

Frond apex 

Widest portion of lamina 
Pinna attachment 
Pinna shape 

Pinna apex 

Sterile v. fertile frond 


Pinnules of fertile fronds 
Indusium length 

Indusium margin 4 ' 

Sporangia 1 stalks 

Spore color 

Spore surface 

Mean spore dimensions 


A . angustum 


erect or ascending, more condensed 
linear lanceolate, 

8-10 X 1*5-2 mm, 
brown to dark brown 
up to half the lamina length 


elliptic or rhombic 
gradually tapered to an acute 
to obtuse base 
acute to acuminate 
near or just below middle 
short-stalked or sessile 


c 


)blong-lanceolate, usually widest 
at middle and not parallel-sided 
acute to acuminate 

tending toward dimorphism, the segments 
of fertile fronds narrower and more 
acute than those of the sterile fronds 
narrowly lanceolate and acute 
tending to he shorter than in asplanioides , 
up to 1.1 mm 

irregularly dentate, and/or eiliate with 
eglandular hairs 

bearing glandular hairs or less often 
secondary sporangia 
yellowish 
sparsely papillate 
38.6 X 24.7 u 


A. asplenioides 

ascending to creeping, more extended 
lanceolate, 3—9 X 2—3 mm, bronze to light 
brown or brown 

equaling lamina length 

narrowly deltoid lanceolate, ovate-lanceolate to lanceolate 
slightly reduced and truncate at base 

acute, acuminate, or more-or-less caudate 
second pinna pair 
usually stalked 

oblong-lanceolate to lanceolate, nearly parallel-sided 
acute 

not tending toward dimorphism 


oblong or linear-oblong and obtuse 

tending to be longer than in angustum, up to 1.3 mm 

eiliate with glandular or non-glandular hairs 
as long as indusial width 
consistently bearing glandular hairs 


brownish-yellow to dark-brown or black 

wrinkled or reticulated exospore, sometimes nigrescent 

36.0 X 25.5 g 


Descriptions of the indusium margin in the two taxa are inconsistent in the literature. The table entry is a tentative consensus. 
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Table 2. Collection localities tor Athyrium filix-femina populations examined in this study. Localities are ordered from north to south. 


Population 

Designation 

Collector(s) and 
collection number 

Year 

Collected 

Locality 

•r 

No. of isozyme 
samples 

No. of spore 
samples 

A. filh 

c-Jemina var. ungustum 





Mt. Sainte Hilaire 

C. Wurth s. n. 

1997 

Province du Quebec, along hiking 
trail through forest. 

Mont Ste.-Hilaire 

26 

7 

Mast Landing 

L. B. Kass s. n. 

1987 

Maine, Cumberland County, Mast 
Landing Sanctuary, Freeport 

0 

4 

North Hudson 

C. Worth and 

C. Caplen s. n. 

1997 

New York, Essex County, Abies/Thuja 
woods and stream 
bank, 11.3 mi VV of North Hudson 
and 1-87, exit 387 
to Blue Ridge 

26 

10 

Schenevus 

C. Worth and 

C. Caplen s. n. 

1997 

New York, Otsego County, woods 
near edge E side 1-88, 
near Schenevus 

11 

1 

Barnet 

C. Werth , D. Conant, 
and C. Caplen s. n. 

1997 

Vermont, Caledonia County, ditches 

m/ 

and woods edge along dirt 
road in vicinity of Conant farm, 

mt 

near Barnet 

72 

20 

Ithaca 

CL Kelloff 560 

1988 

New York, Tompkins County, fen 
along Rt. 96 south of junction 
with Route 13 approximately 2 
miles south of Ithaca 

0 

20 

Binghamton 

CL Kelloff 559 

J m 

1988 

New York, Broome County, floodplain 
of Chenango R., 

1.5 miles N of Binghamton 

0 

20 

Ralston-1 

CL Kelloff 569 

1988 & 

1995 

Pennsylvania, Lycoming County, 
flood plain along Route 

14, 7 miles north of Ralston 

20 

58 

Ralston-2 1 

CL Kelloff 570 

9* m 

1988 

Pennsylvania, Lycoming County, 
flood plain along Route 

14, 3 miles south of Ralston 

0* 

78 


c 

o 

to 
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Table 2. Continued. 


tx 


Population 

Designation 


Collector(s) and 
collection number 


Shirley 


A. fihx-jeminu var. asplenioides 

CL Kelloff 128 o 


Breathed Mountain 


George Mason 


Hopewell 


Mount joy St or 


Pond Drain 


Pipevine Ilollow 


Sandy Run Swamp 


CL Kelloff 25 


CL Kelloff 101 


CL Kelloff 505 


CL Kelloff 564 


C . Werth, 

E. Potter, and 
K. Sc in ret to s. n. 


CL Kelloff s. n. 


CL Kelloff 1205 


Year 

Collected 


1997 


1987 


1988 


1988 & 
1995 


1988 & 
1995 


1997 & 
1998 


1987 


1997 


Locality 


New Jersey, Elmer County, floodplain 
along County Road 
611, 1.5 miles east of jet. NJ State 
Route 77, in Shirley 
West Virginia, Tucker County, 
Breathed Mountains 
Virginia, Fairfax County, creekside on 
George Mason University campus 
in Fairfax 

Virginia, Fauquier County, floodplain 
along Route 629, 3.3 miles N 
of Route 601 near Hopewell 
Virginia, Stafford County, hog at 
Mount joy Store, junction 
of routes 611 and 633 
Virginia, Giles County, north-facing 
wooded slope along rd. 
to White Pine Lodge and above 
Pond Drain, NW of 
rt. 613 and Mountain Lake 
'Tennessee, Davidson County, 

Pipevine Hollow 

North Carolina, Onslow Countv, creek 
bank and roadside ditch E side of 
Sandy Run Swamp along Haws 
Run Road, just W of Padgett, 


( 


:a. 60 air km NL of Wilmington 


No. of isozyme 
samples 


40 


0 


0 


30 


30 


50 


0 


32 


No. of spore 
samples 


27 


6 


30 


42 


63 


0 


?5 


18 


Initial isozyme data were reported for 78 individuals from Ralston-2 (Kelloff 1990), but this population was not resampled during 1995 and iso¬ 
zyme data are not included in the present report (see text for further explanation). 
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Of these, seven were examined only for spores, while for ten (five of each 
taxon), population samples for isozyme electrophoresis were obtained as well. 
At most localities, plants were collected in a transect and individuals num¬ 
bered successively. Voucher specimens for localities collected by CLK are 
deposited at the George Mason University Herbarium (GMUF), those collected 
by CRVV at the E. L. Reed Herbarium of Texas Tech University (TTC). 

Spore Morphology. —Spores were collected from fresh samples gathered from 
16 of the localities listed in Table 2. The spores were examined for surface 
characteristics of the perispore using light microscopy (LM; 409 individuals), 
and scanning electron microscopy (SEM; one to two individuals from ten 
populations, 13 from the Shirley, NJ population). For exine ultrastructure, 
transmission electron microscopy (TEM) was used to examine four individ¬ 
uals, 2 of each taxon. For LM, spores were mounted by gently tapping a 
pinnule bearing mature sporangia over a glass slide and then adding a few 
drops of Permount ,M mounting medium and a glass coverslip. Spores were 
examined under a standard brightfield optic system. 

For SEM, air dried spores were affixed to glass coverslips with ethanol, 
mounted on SEM stubs, and sputter coated with gold or carbon/gold palla¬ 
dium using the Hummer VII sputtering system. The spores were then exam¬ 
ined using either a Hitachi S-530 or a LEO 440 scanning electron 
microscope. From each SEM sample at least 100 spores were examined. 

For TEM, each spore sample was embedded in 1% agar as a pellet and 
fixed in sodium veronal-acetate buffered potassium permanganate (Dawes, 
1971). The sample pellet, cut up into 1.0 mm cubes, was dehydrated in a 
graded series of ethanol. After dehydration, the spores were infiltrated with 
a series of propylene oxide and Spurrs plastic embedding medium in a 
labeled BEEM capsule and polymerized in an oven at 70 C for 16 hours. The 
blocks were then thin-sectioned with a Dupont diamond knife, the sections 
picked up on #200 mesh grids and examined with a JEOL 100G transmisson 
electron microscope. Post-staining of the spore sections was not necessary 
when fixed with the buffered potassium permanganate (Dawes, 1971). 

Isozyme Electrophoresis. —Fronds to be analyzed electrophoretically were 

collected in the field and kept refrigerated in plastic bags until homogenized 
within a few days. For each sample, approximately 1 g leaf tissue was 
ground using eight drops of '■microbuffer" (VVerth, 1985) enhanced with 5% 
PVP-40T and 1% 2-mercaptoethanol. The tissue was ground in a porcelain 
spot plate using a small test tube as a pestle and sand to facilitate grinding. 
The extract was either absorbed onto filter paper wicks and used immedi¬ 
ately or stored at — 78 G for up to four months before being thawed and 
absorbed onto wicks. Starch gel concentration varied from 11% to 13.5% 
depending on the properties of individual starch lots. 

The following eleven enzymes were analyzed using the buffer systems in¬ 
dicated: on lithium hydroxide (Selander et al., 1971)—glutamate oxaloace- 
tate transaminase (GOT), hexokinase (HK), and leucine aminopeptidase 
(LAP); on system number 6 (Soltis et al., 1983)—phosphoglucose isomerase 
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(PGI), phosphoglucomutase (PGM), and triose-phosphate isomerase (TPI); on 
morpholine-citrate pH 8.2 (Werth, 1991)—aldolase (ALD), isocitrate dehydro¬ 
genase (IDH), malate dehydrogenase (MDH), 6-phosphogluconate dehydro¬ 
genase (6PGDH), and shikimate dehydrogenase (SKDH). Gels were run until 
the bromophenol blue marker migrated 10 cm for lithium-hydroxide and 
number 6, or 12 cm for morpholine-citrate. Staining followed Werth (1985) 
as modified for the zvmecicle technique (Werth, 1990). Isozyme banding pat¬ 
terns were readily interpreted as allelic variation for single gene loci from 
whole leaf extracts (cf. Gastony and Darrow, 1983) using standard principles 
(Wendel and Weeden, 1989; Murphy et al„ 1990). Data were analyzed using 
the BIOSYS-1 computer program (Swofford and Selander, 1981), entered into 
the program as individual multilocus genotypes. Occasionally, sets of two 
or three successive samples from the field possessed identical genotypes. As 
Athyrium filix-feminci is known to occasionally form clones up to nearly 
20 m in extent (Sciaretta et al., ms), such sets were assumed to comprise sin¬ 
gle genets and each set was represented by only a single entry in the data set. 

Initial isozyme data were obtained from four populations, two of each 
taxon (data reported in Kelloff, 1990). During 1995, three of these popula¬ 
tions (Ralston-1, PA, Hopewell, VA, and Mountjoy Store. VA) were revisited 
and re-sampled, and an expanded isozyme data set obtained; this was aug¬ 
mented by sampling of additional populations during 1997 and 1998. The 
initial data were largely consistent with the more recent results; however, 
only the latter are reported here to avoid including inadvertently re-sampled 
individuals. 


Results 

Spores, Light Microscopy. —Observation of spores through the light micro¬ 
scope (LM) revealed features largely consistent with previous reports of dif¬ 
ferences between the two taxa (Butters, 1917; Lugardon, 1971; Schneller, 
1989). Spores of both A. angustum (Fig. 2a, b) and A. asplenioides (Fig. 3a. b) 
were monolete with an unbranched longitudinal proximal laesura, ovoid to 
ellipsoidal in polar view, and possessed a perispore (terminology of Erdtman. 
1969). Although most of the surface details were obscure when viewed under 

LM, the difference between the two taxa was readilv seen in their outline. The 

•/ 

spores of A. angustum had a somewhat smooth surface (Fig. 2a, b), while those 
of A. asplenioides had an uneven surface (Fig. 3a, b). 

Spore color has been described as differing between A. angustum with vel- 
lowish spores, and A. asplenioides with yellowish-brown to blackish spores 
(Table 1). This characterization evidently is based on appearance of mass 
spores under reflected light as viewed unaided or with a dissecting micro¬ 
scope. When viewed under LM using transmitted light, spores of both taxa 
appeared yellowish, with the exception of occasional A. asplenioides spores 
that displayed an especially dark or blackish reticulation overlaving the yel¬ 
low surface of the spore. 
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Mt. Sainte Hilaire 

-Barnet 
-North Hudson 

-Schenevus 
Ithaca 

Binghamton 

-Ralston-1 

-Ralston-2 

Shirley 

Breath Mountain 
Hopewell 
George Mason 
Mountjoy Store 
Pond Drain 
Pipevine Hollow 
Sandy Run Swamp 


A. angustum 
A. asplenioides 



500 Miles 



Hr,. 1. Map showing overlapping ranges ot two Athynrium taxa, angustum and asplenioides, in 
eastern North America and location of populations studied. 


Scanning Electron Microsc:opy. —Observations of spores under SEM were 
consistent with those ot Tryon and Lugardon (1991) and Schneller (1989) in 
revealing that the perispore ot A. angustum (Fig. 2c—f) differs markedly from 
that of A. asplenioides (Fig. 3c—f). The perispore of A. angustum was papil¬ 
lose. its surface completely and evenly covered with minute verruciform 
(i.e., wart-like) projections. On some spores of A. angustum, thin, irregular 
laminae of perispore material appeared as occasional “flakes.” In contrast, 
the perispore ot A. asplenioides was rugose, characterized by a network of 
nniri (i.e. ridges/walls) that meandered and joined to surround lacunae. The 
surface ot the lacunae completely lacked the projections as in A. angustum. 

Transmission Electron Microscopy.— Under TEM the similarities and differ¬ 
ences between the structural and sculptural features of the spore wall of the 

y l ^ A. angustum and A. asplenioides (Figs. 

2t and 3f) possessed a thick external exospore (Ee). Although not evident in 
these sections, the tern exospore is composed of feuillettes, i.e., fused plates 
(Lugardon, 1976, 1979; Tryon, 1986). TEM revealed that the sculptural ele¬ 
ments composing the surface details seen under SEM were not derived from 
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Fk;. 2. Spores of Athyrium angustum . A, B. Spores viewed under LM. Note smooth appearance 
of perispore. C. D, E. Spores viewed under SEM. Note papillate perispore surface. F. Spore wall 
viewed under TEM. Note low perispore (p), external exospore (ee), and internal exospore (ei). 


the exospore, as for example in the spores of Trichomcines or Protomarattia 
(Tryon, 1986), but instead from the perispore, the outer stratum circumjacent 
to the external exospore. The perispore of both taxa fit tightly to the external 
exospore. This stratum in the spores of A. angustum (Fig. 2f) was a shallow. 
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Fig. 3. Spores of Athyrium asplenioides. A, B. Spores viewed under LM. Note rugose, reticulate 

appearance of perispore. C, D, E. Spores viewed under SEM. Note reticulating inflated folds of 

perispore surface. F. Spore wall viewed under TEM. Note inflated fold in perispore (p). external 
exospore (ee), and faint internal exospore (ei). 


fairl\ uniform layer. In A . asplenioides spores, what appeared to be muri 
under the SEM (Figs. 3c-e) were more clearly seen under TEM to be inflated 
folds in the perispore because the columellae, characteristic of true muri, 
were lacking (Fig. 3f). Beneath the external exospore on A . ongustum was 
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a thin endospore or internal exospore (Ei), a layer that is deposited shortly 
before germination of the spore (Tryon, 1986). A remnant of this layer was 
also visible in A. asplenioides. 


Taxon Fidelity of Spore Morphotypes and Evidence for Hybridization. —All 
populations of A. angustum sampled for this study exhibited the smooth, 
micropapillate spore morphotype, and all populations of A. asplenioides 
except one exhibited the rough, rugose spore morphotype. Reliance on light 
microscopy, which may fail to reveal the subtle detail of the perispore, and 
on color, which varies in A. asplenioides despite the overall similarity of 
perispore in this taxon, may have led to some confusion in the previous lit¬ 
erature. For example, Liew (1971) concluded that most individuals of both 
A. angustum and A. asplenioides produced both kinds of spores, entirely in¬ 
consistent with our observations. 

The exceptional A. asplenioides population was the northernmost popula¬ 
tion sampled, Shirley, NJ. At this site the spores had perispore sculpturing 
characteristics of both A. angustum and A. asplenioides. While none of the 
plants from this site possessed A. angustum- type spores, spore arrays of the 
13 plants examined under SEM possessed unique and variable morphologies 
that suggested an influence from A. angustum (Fig. 4). The greatest degree of 
A. angustum influence was exhibited by six individuals in which prevailing 
perispore morphotype was very irregular in appearance, combining an array 
of sculptural elements that were difficult to characterize (Fig. 4a, b). These 
included irregularly shaped papillae, briefly-extending isolated rnuri, and a 
preponderance of irregular small lamellae reminiscent of the "flakes - ’ seen 
in A. angustum. This prevailing morphotype is readily interpreted as inter¬ 
mediate between the highly distinct morphotypes of A. angustum and A. as¬ 
plenioides, seeming to average their disparate surface features and indicating 
that these individuals may be hybrids between the two taxa. Many of the 
highly variable spores of these putative Fl hybrids were not entirely inter¬ 
mediate, but rather showed tendencies toward either the smooth A. angus¬ 
tum type or the rough A. asplenioides type (Fig. 4a, c-g). Spores that fully 
resembled parental types were discovered in the spore arrays of these six in¬ 
dividuals, although these extremes were decidedly rare. The spores of these 
putative hybrids were not abortive, as is typically the case in interspecific 
fern hybrids, but rather appeared normal, i.e., well-filled, and were found to 
be viable as evidenced by successful germination (unpublished data). More¬ 
over, the seven other individuals in the Shirley. Nf, population examined 
under SEM possessed spore arrays exhibiting various degrees of intermediacv 
but strongly skewed toward the A. asplenioides morphology. These are 
hypothesized to include backcrosses between first-generation hybrids and 
A. asplenioides (two individuals) and later generation backcrosses (five indi¬ 
viduals). 


Allozymes, General. —The eleven enzymes assaved were coded bv 17 inter- 

v i/ J 

pretable loci of which only one ( Pgi-1 ) was invariant across all ten popula¬ 
tions. The remaining sixteen loci [Aid, Got, Hk. Idh-1, Lap, Mdh-1, Mdh-2, 
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Fir,. 4. Spores of putative hybrids between Athvrium angustum and A. asplenioides viewed 

l _ » * i, # A* 

under SEM. A—G, spores ot putative first-generation hybrid. A. Array of spores. Note variable 
morphology and preponderance of “flaky” intermediates. B. Close-up of “flaky” intermediate 
spore morphotype. C—G. Examples ot spore morphotvpes occurring in putative F-l hybrid, vary¬ 
ing from A. angustu rn-l ike (C) to A. asplenioides-\ike (G). H—M, spores of putative backcross 
between F-l hybrid and A. asplenioides. H. Array of spore morphotvpes exhibited bv the puta¬ 
tive backcross. Note prevalent of A. asplenioides -like spores. I. Close-up of A. asplenioides -like 
spore from backcross. Note irresolute nature of inflated folds and tendency toward "flakiness”. 

J-M. Examples of spore morphotvpes occurring in putative backcross. Note variation from 
"flaky” intermediate (J) to A. asplenioides -like (M). 
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Mdh-3, Mdh-4, Pgi-2, Pgm-2, 6-Pgd-l, 6Pgd-2, Skdh, Tpi-1, and Tpi-2) were 
variable in at least one population, and the frequencies of the alleles were 
computed (Table 3). There was a strong tendency for all populations to share 
the more common alleles and some of the less frequent ones as well. Most of 
the loci were weakly polymorphic, a single allele predominating in all popu¬ 
lations of both A. angustum and A. asplenioides. In contrast, at the four 
most polymorphic loci ( Idh-1, Pgi-2, Pgm-2, and Tpi-2), population allele fre¬ 
quencies were similar within A. angustum and A. asplenioides respectively, 
but strikingly different between the two taxa (Table 3). 

The most contrasting locus was Idh-1 for which alleles Idh-1'* and Idh-1 1 


prevailed in A. angustum at frequencies ranging from 0.545 to 0.676 and 
0.296 to 0.455 respectively, while allele Idh-1 A was the most frequent in 
A. asplenioides with frequencies ranging from 0.828 to 0.950. At Pgi-2 both 
taxa shared Pgi-2 as their most frequent allele; however, Pgi-2 was repre¬ 
sented in all A. angustum populations at much higher frequencies, ranging 
from 0.154 to 0.450. than in populations of A. asplenioides for which the fre¬ 
quency of this allele ranged from 0.000 to 0.106. At Pgm-2, both taxa shared 

n 

Pgm-2’ and Pgm-2' as principal alleles, the former prevailing in A. angus¬ 
tum populations at frequencies ranging from 0.856 to 1.000. the latter pre¬ 
vailing in A. asplenioides populations at frequencies ranging from 0.580 to 
0.857. At Tpi-2 both taxa shared three prevalent alleles, with Tpi-2 A being of 
substantial frequency in all populations. Allele Tpi-2 h occurred at high fre¬ 
quencies, ranging from 0.355 to 0. 550, in A. angustum populations whereas 

mm*. * . i ' « 

Tpi-2 ' was infrequent in this taxon, occurring at frequencies from 0.000 to 

_ f » 

0.100. Conversely, Tpi-2 often was the most frequent allele in A. asple¬ 
nioides populations, occurring at frequencies from 0.316 to 0.606, while Tpi- 

fl _ _ 1 

2 was rarer, occurring at frequencies from 0.000 to 0.255. The contrasting 
allele frequency trends for these four loci were highly consistent among pop¬ 
ulations within each taxon. Exceptions to these frequencies trends were in 
the northernmost sampled A. asplenioides population. Shirlev, from south¬ 
ern New Jersey for the characteristically A. angustum alleles Idh-1 B and Pg/- 
2 ’ and likewise for Tpi-2’ in the highest elevation (1300 m) population of 
A. asplenioides sampled. Pond Drain, from the mountains of southwestern 
Virginia. Moreover, five of the six Shirley, NJ. individuals hypothesized to 
be first-generation hybrids on the basis of spore morphology were hetero¬ 
zygous for A. angustum and A. asplenioides marker alleles for at least three 
of the four most divergent loci. i.e.. Idh-1 AH , Pgi-2 h( \ Pgm-2 B( (scored for 
only two of these individuals) and Tpi-2 . No other individuals in the entire 
data set possessed this genotype combination. 


Genetic Variation. —Genetic variation was quantified for each population 
and the species as a whole by computing three standardly used indices. Val¬ 
ues for percent loci polymorphic (P) ranged from 23.5% to 47.1%, with a 
mean of 36.48%r, for mean number of alleles per locus (A), the range was 1.5 
to 2.5, mean 1.97; and for mean expected heterozygosity (He), the range was 
0.112 to 0.147, mean 0.129 (Table 4). The mean values for these indices in 




Table 

; 3. Allele 

frequencies 

for 17 isozyme 

loci in ten 

populations of Athyrium filix-femina 

in eastern 

North America. 


Locus 

Allele 



angustum 





cisplenioides 


Mt. Ste. 

1 Ida ire QL 

Barnet 

VT 

North 

Hudson NY 

Scheneviis 

NY 

Ralston 

PA 

Shirley 
NJ ' 

Mount joy Hopewell 
Store VA VA 

Sandy Run Pond Drain 

mr 

Swamp NC VA 

Aid 

A 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.969 


B 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.021 


C 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.010 


(N) 

19 

29 

26 

11 

1 

8 

1 

1 

32 

48 

Got 

A 

0.000 

0.033 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


B 

1.000 

0.951 

1.000 

1.000 

1.000 

0.955 

1.000 

0.917 

1.000 

0.980 


c 

0.000 

0.016 

0.000 

0.000 

0.000 

0.045 

0.000 

0.083 

0.000 

0.010 


D 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.010 


(N) 

10 

61 

5 

6 

20 

33 

30 

30 

12 

50 

Hk* 

A 

0.000 

0.000 

0.000 

0.000 

0.150 

0.000 

0.017 

0.000 

0.000 

0.016 


B 

0.960 

0.965 

1.000 

0.955 

0.800 

0.868 

0.967 

0.967 

0.893 

0.952 


c 

0.040 

0.023 

0.000 

0.045 

0.050 

0.105 

0.017 

0.017 

0.089 

0.032 


E 

0.000 

0.012 

0.000 

0.000 

0.000 

0.025 

0.000 

0.017 

0.000 

0.000 


F 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.018 

0.000 


(N) 

25 

43 

19 

11 

10 

19 

30 

30 

28 

31 

ldh-1* 

A 

0.019 

0.007 

0.000 

0.000 

0.000 

0.828 

0.950 

0.833 

0.859 

0.920 


B 

0.596 

0.676 

0.596 

0.545 

0.650 

0.172 

0.017 

0.133 

0.094 

0.060 


C 

0.385 

0.296 

0.385 

0.455 

0.350 

0.000 

0.000 

0.000 

0.000 

0.020 


I) 

0.000 

0.021 

0.019 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


F 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.033 

0.033 

0.047 

0.000 


(N) 

26 

71 

26 

11 

20 

32 

30 

30 

32 

50 

La p 

A 

0.000 

0.000 

0.000 

- 

0.000 

0.000 

0.018 

0.000 

0.000 

0.000 


B 

0.100 

0.038 

0.063 

- 

0.175 

0.000 

0.036 

0.000 

0.031 

0.016 


C 

0.850 

0.885 

0.750 

• 

0.550 

1.000 

0.946 

1.000 

0.938 

0.935 


D 

0.050 

0.058 

0.188 

- 

0.275 

0.000 

0.000 

0.000 

0.016 

0.032 


E 

0.000 

0.019 

0.000 

• 

0.000 

0.000 

0.000 

0.000 

0.016 

0.000 


F 

0.000 

0.000 

0.000 

- 

0.000 

0.000 

0.000 

0.000 

0.000 

0.016 


(N) 

10 

26 

16 

0 

20 

33 

28 

23 

32 

31 
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Tahle 


angustuni 


Locus 

Allele 

ML Ste. 

1 lilaire QE 

Barnet 

VT 

North 

Hudson NY 

Schenuvi 

NY 

Mdh-l 

A 

0.981 

0.950 

0.942 

1.000 


B 

0.019 

0,020 

0.058 

0.000 


c 

0.000 

0.030 

0.000 

0.000 


(N) 

26 

50 

26 

11 

Mdh-2 

A 

1.000 

0.990 

1.000 

1.000 


B 

0.000 

0.010 

0.000 

0.000 


(N) 

26 

50 

26 

11 

Mdh-3 

A 

1.000 

1.000 

1.000 

1.000 


B 

0.000 

0.000 

0.000 

0.000 


(N) 

26 

50 

26 

11 

Mdh-4 

A 

0.960 

0.920 

0.962 

1.000 


B 

0.040 

0.070 

0.038 

0.000 


c 

0.000 

0.010 

0.000 

0.000 


(N) 

25 

50 

26 

11 

Pgi-1 

A 

1.000 

1.000 

1.000 

1.000 


(N) 


67 

21 

11 

Pgi-2 

A 

0.000 

0.000 

0.000 

0.000 


B 

0.000 

0.014 

0.019 

0.000 


C 

0.019 

0.007 

0.000 

0.000 


D 

0.000 

0.000 

0.000 

0.000 


E 

0.HM2 

0.715 

0.808 

0.500 


F 

0.019 

0.014 

0.000 

0.000 


G 

0.269 

0.250 

0.154 

0.409 


H 

0.000 

0.000 

0.000 

0.000 


i 

0.000 

0.000 

0.019 

0.000 


i 

0.000 

0.000 

0.000 

0.091 


(N) 

26 

72 

26 

11 

Pgm-2 

A 

0.024 

0.000 

0.000 

0.000 


B 

0.857 

0.856 

0.906 

1.000 


c 

0.119 

0.122 

0.094 

0.000 


t. Continued. 


cisplenioides 


Ralston 

PA 


Shirley 
N| * 


Mountjoy I lopewe 
Store VA VA 


Sandy Run Pond Drain 


Swamp NC 


VA 


0.967 


1.0( H ) 

0.000 

0.000 

20 

1.000 

0.000 

20 

1.000 

0.000 

20 

1.000 

0.000 

0.000 

13 

1.000 

to 

0.000 

0.000 

0.000 

0.000 

0.550 

0.000 

0.450 

0.000 

0.000 

0.000 

20 

0.000 

0.975 

0.025 


0.939 

0.061 

0.000 

33 

1.000 

0.000 

33 

0.985 

0.015 

33 

1.000 

0.000 

0.000 

33 

1.000 

13 

0.015 

0.000 

0.000 

0.045 

0.742 

0.045 

0.106 

0.045 

0.000 

0.000 

33 

0.000 

0.143 

0.857 


0.033 

0.000 

15 

0.933 

0.067 

15 

0.967 

0.033 

15 

0.967 

0.033 

0.000 

15 

1.000 

25 

0.033 

0.017 

0.000 

0.100 

0.800 

0.017 

0.033 

0.000 

0.000 

0.000 

30 

0.000 

0.300 

0.683 


1.000 

0.000 

0.000 

20 

0.950 

0.050 

20 

1.000 

0.000 

20 

0.950 

0.050 

0.000 

20 

1.000 

24 

0.000 

0.133 

0.017 

0.000 

0.850 

0.000 

0.000 

0.000 

0.000 

0.000 

30 

0.083 

0.100 

0.800 


0.969 

0.000 

0.031 

32 

1.000 

0.000 

32 

0.953 

0.047 

32 

1.000 

0.000 

0.000 

32 

1.000 

27 

0.000 

0.000 

0.000 

0.078 

0.906 

0.000 

0.016 

0.000 

0.000 

0.000 

32 

0.000 

0.222 

0.704 


1.000 

0.000 

0.000 

50 

1.000 

0.000 

50 

1.000 

0.000 

50 

1.000 

0.000 

0.000 

50 

1.000 

38 

0.000 

0.000 

0.010 

0.010 

0.930 

0.020 

0.030 

0.000 

0.000 

0.000 

50 

0.020 

0.390 

0.580 





H 

a 

2! 

Z 






















1 ~ 


Tablk 3. Continued. 


angustum cisplenioides 


Locus 

Allele 

Mt. Ste. 
Hilaire QE 

Barnet 

VT 

North 

Hudson NY 

Schenevus 

NY 

Ralston 

PA 

Shirley 

NJ ' 

Mountjoy 
Store VA 

Hopewell 

VA 

Sandy Bun 

Swamp NC 

Pond Drain 

VA 



1 ) 

0.000 

0.022 

0.000 

0.000 

0.000 

0.000 

0.017 

0.017 

0.074 

0.000 



E 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.010 



(N) 

21 

45 

16 

11 

20 

21 

30 

30 

27 

50 


6Pgd-1 

A 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 



B 

1.000 

1.000 

1.000 

1.000 

1.000 

0.957 

1.000 

1.000 

1.000 

1.000 



c 

0.000 

0.000 

0.000 

0.000 

0.000 

0.043 

0.000 

0.000 

0.000 

0.000 

s' 


(N) 

11 

38 

11 

6 

10 

23 

14 

15 

15 

44 

3 

6Pgd-2 

A 

0.000 

0.000 

0.000 

0.000 

0.000 

0.015 

0.033 

0.017 

0.017 

0.021 

L J 


B 

1.000 

1.000 

1.000 

1.000 

1.000 

0.985 

0.950 

0.967 

0.983 

0.947 

*7 


C 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.017 

0.017 

0.000 

0.021 



D 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.011 

H 


(N) 

26 

69 

21 

11 

20 

33 

30 

30 

30 

47 

2: 

Skdh 

A 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.067 

0.000 

0.000 

c 


B 

0.981 

0.955 

0.942 

0.955 

1.000 

1.000 

1.000 

0.933 

1.000 

0.990 

$3 


C 

0.019 

0.036 

0.058 

0.045 

0.000 

0.000 

0.000 

0.000 

0.000 

0.010 

2: 

*> 


i) 

0.000 

0.009 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

• • 


(N) 

26 

56 

26 

11 

20 

33 

30 

30 

32 

50 


Tpi - 7 

A 

0.947 

1.000 

0.981 

0.900 

1.000 

0.970 

1.000 

1.000 

0.984 

1.000 

u 

L 


B 

0.053 

0.000 

0.000 

0.100 

0.000 

0.030 

0.000 

0.000 

0.016 

0.000 

v —1 


c 

0.000 

0.000 

0.019 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 



(N) 

19 

47 

26 

5 

10 

33 

25 

25 

32 

18 

is: 

Tpi-2* 

A 

0.579 

0.601 

0.404 

0.400 

0.450 

0.318 

0.440 

0.420 

0.406 

0.429 



B 

0.368 

0.355 

0.538 

0.500 

0.550 

0.076 

0.100 

0.000 

0.063 

0.255 



C 

0.053 

0.036 

0.058 

0.100 

0.000 

0.606 

0.460 

0.580 

0.531 

0.316 

5 


E 

0.000 

0.007 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

50 


(N) 

19 

69 

26 

5 

10 

33 

25 

25 

32 

49 



is: 


ro 
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these eastern North American Athvrium populations were somewhat greater 
than the means tor ferns (P=36.0%, A = 1.65, H=0.109) obtained by averaging 
across 32 taxa (Li and Haufler, 1999), which are in turn similar to the means 
for angiosperms (P = 34.2%, A = 1.53, H E = 0.113; Hamrick and Godt. 1989). 
Mean values across A. angustum populations for P (34.12%) and A (1.84) 
were slightly lower than those for A. asplenioides (mean P=38,84%, mean 
A = 2.1), while the mean value in A. angustum for H E =0.138 was greater than 
that of 0.119 in A. asplenioides. This indicates that while A. asplenioides 
possesses greater allelic diversity, frequencies of alleles within loci are more 
evenly distributed in A. angustum populations. 

Values tor mean observed heterozygosity (Ho) were also computed for 
comparison to H E . Values of H () were very similar to (although tending to be 
slightly greater than) those for H E , suggesting a tendency toward random 
mating. 


Comparison to Hardy-Weinberg. —Genotype proportions for each polymorphic 
locus in each population were compared to Hardy-Weinberg expected values 
by computing the fixation index F, and the statistical difference of F from 0 
was evaluated using the chi-square test, with pooled genotype classes if the 
number of alleles exceeded 2 ( fable 5). The test was considered valid if 
two of the three genotype classes were represented by expected values 5= 5. 
Of 42 validly tested loci, 40 conformed to Hardy-Weinberg expectations. 
Moreover. 45 out of 51 non-valid tests also indicated conformance to Hardv- 
Weinberg values, despite the tendency for non-valid tests to indicate false 
non-conformance due to low expected values. Thus, mating was inferred to 
approximate random mating via predominant outcrossing between gameto- 
phytes. as appears to be the general case in most ferns (Soltis et al ., 1988) 
including Athyrium species (Schneller, 1979). 


Genetic Relatedness of Populations and Taxa. —The degree of genetic diver¬ 
gence among populations ( fable 6) was quantified by computing F-statistics 
(Wright, 1965. 1978), including hierarchical analysis (Wright. 1978). Values 
for Fst computed among all populations varied among loci from a low of 
0.019 for 6Pgd-2 to high values of 0.468 for Pgm-2 and 0.460 for ldh-1, 
the latter two loci having exhibited the greatest allele frequency difference 
between A. angustum and A. asplenioides. The value across loci of F. ST = 
0.255 indicated very substantial differentiation among populations. This val¬ 
ue is very high in comparison to other fern species examined, exceeding, for 
example, computed values for mean F S t of 0.024 among populations of Polv- 
stichum munitum ranging from Oregon to Idaho (Soltis et al.. 1987), F ST = 
0.152 among Pteridium aquilinum populations ranging from Massachusetts 
to Florida (Speer et al., 1998). and F S t = 0.100 to 0.248 in various species of 
Dmopteris ranging widely across eastern North America (Werth. ms.). 

To evaluate the contribution of differences between A. angustum and A. 
asplenioides to overall population differentiation, hierarchical F-statistic 
analysis (Wright, 1978) was carried out (Table 6). For most individual loci, 
as well as for the combined values across loci, the variance between the two 




I able 4. Estimates oi genetic variation at 17 loci 

in ten populations of 

Athyrium filix-femina s. 

/. (standard errors in 

parentheses). 

Population 

Mean sample size 
per locus 

Mean no. of alleles 
per locus (A) 

Percent loci 
polymorphic (P) 

Mean heterozygosity 

Hardy-Weinberg 

Observed (HJ Expected (He) 

A . august uni 

Mt. Ste. Hilaire, QE 

21.4 (1.4) 

1.9 (0.2) 

35.3 

0.150 (0.056) 

0.140 (0.047) 

Barnet, VT 

52.5 (3.4) 

2.5 (0.3) 

41.2 

0.134 (0.042) 

0.141 (0.042) 

North Hudson, NY 

21.7 (1.5) 

1.8 (0.2) 

41.2 

0.133 (0.046) 

0.136 (0.047) 

Schenevus, NY 

9.1 (0.8) 

1.5 (0.2) 

23.5 

0.145 (0.071) 

0.126 (0.055) 

Ralston, PA 

15.5 (1.4) 

1.5 (0.2) 

29.4 

0.150 (0.061) 

0.147 (0.057) 

Mean 

24.04 

1.84 

34.12 

0.142 

0.138 

A. asplanioides 

Shirley, NJ 

28.2 (2.0) 

2.0 (0.3) 

35.3 

0.127 (0.041) 

0.127 (0.040) 

Mount joy Store, VA 

22.8 (2.1) 

2.2 (0.3) 

41.2 

0.133 (0.052) 

0.123 (0.042) 

Hopewell, VA 

23.7 (1.8) 

1.9 (0.2) 

47.1 

0.125 (0.038) 

0.118 (0.036) 

Sandy Run Swamp, NC 

28.9 (1.5) 

2.0 (0.2) 

35.3 

0.104 (0.039) 

0.117 (0.041) 

Pond Drain, VA 

44.5 (2.3) 

2.4 (0.3) 

35.3 

0.118 (0.049) 

0.112 (0.046) 

Mean 

29.62 

2.1 

38.84 

0.121 

0.119 

Mean across all 

26.83 

1.97 

36.48 

0.131 

0.129 

populations 


a A locus is considered polymorphic if the frequency of the most common allele does not exceed 0.95. 
h Unbiased estimate (Nei, 1978). 
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Table 5. Values of the fixation index F as a measure of the conformance of population genotype ratios to those expected under Hardy-Weinberg 
equilibrium. F was computed for each polymorphic locus in each population, and tested for statistical difference from 0 (i.e. conformance to Hardy- 

mm m. am m V a. ^ a m. _ —. _ — __ 


classes. Chi-square 

tests were considered valid only if two of the three genotypic classes were re 

presented by an expected values > 5. Results from 

non-valid tests are 

provided in brackets. Values marked as “ns” 

represent 1 

oci found 1 

o confirm to Hardy-Weinberg ex 

pected proportions (p ^ 0.05), 

and therefore are not statistically 

different from 0. Values marked 

with asterisks are statistically different from 0 at probabilities p ^ 0.05 (one asterisk), 

p ^ 0.01 (two asterisks), or p ^ 

0.001 (three asterisks). 








Mt. Ste. 





Mou nt joy 


Sandy Run 



Hilaire 

Barnet 

N Hudson Schenevus 

Ralston 

Shirley 

•r 

Store 

Hopewell 

Swamp 

Pond Drain 

Locus 

_QK 

VT 

NY NY 

PA 

N] 

VA 

VA 

NC 

VA 

Aid 









(-0.025 ns] 

Got 


-0.040 ns 



(-0.048 ns| 


-0.091 ns 


[-0.015 ns) 

Hk 

1-0.042 ns| 

[—0.028 ns 

[-0.048 ns| 

|0.403 ns) 

-0.124 ns 

1-0.026 ns] 

1-0.026 ns) 

0.266 ns 

[-0.039 nsj 

Idh 

-0.007 ns 

0.102 ns 

0.148 ns (0.267 ns] 

-0.099 ns 

-0.208 ns 

|0.653***| 

-0.163 ns 

0.002 ns 

—0.070 ns 

Lap 

[10.132 ns) 

— 0.087 ns 

0.216 ns 

0.154 ns 


[-0.043 ns| 


10.216*** 

| [-0.046 ns] 

Mdh- 1 

1-0.020 ns) 

0.376* * 

(— 0.061 ns| 


(-0.065 ns| 

[-0,034 ns| 


[-0.032 ns| 


Mdh-2 


(-0.010 ns 




[-0.071 ns] 

-0.053 ns] 



Mdh-3 





(-0.015 ns] 

1-0.034 ns] 


(0.650*** 


Mdh-4 

|1.000***1 

0.192 ns 

(-0.040 ns) 



[-0.034 ns) 

[—0.053 ns] 



Pgi-2 

-0.375* 

0.021 ns 

-0.190 ns 1-0.266 ns| 

-0.212 ns 

-0.054 ns 

0.136 ns 

-0.156 ns 

-0.088 ns 

-0.048 ns 

Pgm-2 

0.240 ns 

0.120 ns 

(-0.103 ns) 

( — 0.026 nsj 

-0.167 ns 

-0.280 ns 

0.028 ns 

0.259 ns 

-0.057 ns 

6Pgd * i 





11 . 000 ***] 





tiPgd-2 





1-0.015 ns) 

1-0.040 ns) 

[-0.026 ns] 

1-0.017 ns] 

-0.038 ns 

Skdh 

[-0.020 ns| 

-0.039 ns 

1-0.061 ns] |-0.048 ns] 




[-0.071 nsj 


1-0.010 ns| 

No. tests 

315] 

H[2| 

4|5[ 0(5] 

3 [ 3] 

5(5| 

2|7] 

5)5| 

5| 31 

5[5| Total: 40[45| 

showing 










conformance 










to HVV 1 










No. of tests 

1(11 

1|0| 

0|0| 0|1| 

0(0) 

0)11 

0|1| 

o(o| 

0|2| 

0(0] Total: 2|6] 

showing 










non-confomuinct 

> 









to HW' 











number o( non-valid lasts in brackets. 
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1 \hi i t>. F statistit analysis (Wright, 1965, 1978), including hierarchical analysis (Wright, 1978), for 16 polymorphic loci across Athvrium filix- 
frmina populations in eastern North America. Statistical difference of F S t from 0 was evaluated using contingency chi-square analysis. For 
hieran Ideal analysis,. I thvvium populations were assigned to their respective taxa (angustum or usphnioides). Differences between values of Fst and 
F local it v t. it 1 1 are attributable to differences in computational method (discussed in S wofford and Selander, 1981), 


Locus 


F Statistic 

s 


Hierarchical F-statistic Analysis (variance 

•r 

components in parentheses) 


Fis 

Fit 

Fst 

F Locality/Total 

l 1 Local ity/Taxon 

F* 

laxon/Total 

Total Limiting 
Variance 

Aid 

-0.025 

-0.002 

0.022 ns 

0.012 (0.00007) 

0.014 (0.00009) 

-0.002 

(-0.00002) 

0.00623 

Got 

-0.060 

-0.017 

0.041 ns 

0.028 (0.00109) 

0.030 (0.00116) 

-0.002 

( 0.00007) 

0.03895 

Hk 

0.116 

0.160 

0.049*** 

0.021 (0.00273) 

0.031 (0.00398) 

-0.010 

(-0.00125) 

0.12816 

Idh-1 

0.043 

0.485 

0.460*** 

0.452 (0.29124) 

0.002 (0.00073) 

0.451 

(0.29050) 

0.64377 

Lap 

0.077 

0.204 

0.137*** 

0.113 (0.02383) 

0.092 (0.01886) 

0.024 

(0.00497) 

0.21012 

Mdh-t 

0.036 

0.060 

0.025 ns 

0.008 (0.00037) 

0.014 (0.00068) 

-0,006 

(-0.00031) 

0.04939 

Mdh-2 

-0.059 

-0.013 

0.044** 

0.017 (0.00042) 

0.016 (0.00040) 

0.001 

(0.00001) 

0.02501 

Mdh-3 

0.302 

0.321 

0.028 ns 

0.007 (0.00013) 

0.002 (0.0003) 

0.005 

(0.00010) 

0.01889 

Mdh-4 

0.208 

0.230 

0.029 ns 

0.010 (0.00046) 

0.014 (0.00067) 

-0.005 

(-0.00021) 

0.04724 

Pgi-2 

-0.143 

-0.012 

0.115*** 

0.096 (0.03914) 

0.040 (0.01551) 

0.058 

(0.02363) 

0.40775 

Pgm-2 

-0.001 

0.469 

0.468*** 

0.459 (0.23433) 

0.045 (0.01293) 

0.434 

(0.22140) 

0.51042 


1.000 

1.000 

0.039 ns 

0.018 (0.00016) 

0.023 (0.00020) 

-0.004 

(-0.00004) 

0.00866 

6Pt>d-2 

-0.032 

-0.013 

0.019 ns 

0.011 (0.00037) 

0.000 (0.00000) 

0.011 

(0.00037) 

0.03325 

Skdh 

-0.052 

-0.019 

0.031*** 

0.013 (0.00061) 

0.011 (0.00054) 

0.001 

(0.00006) 

0.04782 

Tpi-1 

-0.070 

-0.020 

0.047 ns 

-0.003 (-0.00014) 

0.000 (0.00001) 

-0.004 

(-0.00015) 

0.04269 

Tpi-2 

-0.214 

-0.021 

0,159*** 

0.137 (0.08853) 

0.013 (0.00763) 

0.125 

(0.08091) 

0.64843 

Combined across loci 

-0.052 

0.217 

0.255*** 

0.238 (0.68333) 

0.028 (0.06343) 

0.216 

(0.601990) 
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taxa with respect to the total (F XY = 0.216) was an order of magnitude greater 
than variance among populations within taxa (F XY = 0.028). Thus, differences 
between taxa explained most of the variance among populations with respect 
to the total (F XY = 0.238). This result is consistent with and explained by the 
large allele frequency differences at the four most polymorphic loci ( Idh-1. 
Pgm-2, Pgi-2 , and Tpi-2 ) between populations of different taxa as opposed to 
populations of the same taxon. 

Values for Nei's Genetic Identity, 1, (Nei, 1978) and Rogers’ Genetic Simi¬ 
larity, S. (Rogers, 1972) were computed for each pair of populations (Tables 
7 and 8). Values for these indices were consistently higher between popula¬ 
tions of the same taxon, ranging from 1 = 0.990 to 1.000 and S = 0.930 to 
0.975, than between populations of different taxa, ranging from 1 = 0.875 to 
0.938 and S-0.803 to 0.881 (Table 8). 

Populations were clustered using the Unweighted Pair-group Method with 
Averaging (UPGMA) based on both S and I. The two indices resulted in very 
similar dendrograms that differed only in the association among some of the 
A. asplenioides populations; only the dendrogram based on S is illustrated 
(Fig. 5). Two clusters, each comprising all the populations of one taxon, 
were joined at S = 0.849. Within the A. angustum cluster, the two northern¬ 
most populations Mt. Ste. Hilaire, QUE and Barnet, VT were placed as most 
similar, joined at S = 0.975, and to this cluster the other three A. angustum 
populations were joined successively in order from north to south; the 
southernmost A. angustum population Ralston-1, PA joined the A. angustum 
cluster at S = 0.938. The topology of this A. angustum cluster was identical 
in the dendrogram based on I (not shown). In A. asplenioides, the most simi¬ 
lar populations based on S=0.966 were the southernmost population Sandy 
Run Swamp, NG and the next most southeastern occurring population 
Mountjoy Store, VA, and to these were joined the Pond Drain, VA popula¬ 
tion from the mountains of southwestern Virginia at S = 0.958 to form a 
subcluster. A second subcluster, comprising the two northernmost A. asple¬ 
nioides populations Shirley. NJ and Hopewell, VA, joined the more southern 
subcluster at S=0.947. The topology of the A. asplenioides cluster differed 
somewhat in the dendrogram based on I (not shown) indicating that the geo¬ 
graphic ' signal’’ in the A. asplenioides data is weaker than in the A. angus¬ 
tum data. 


Discussion 


The A. filix-femina complex, distributed across four continents and com¬ 
prising as many as four North American taxa with overlapping ranges, pro¬ 
vides an especially suitable context for exploring patterns and processes of 
divergent evolution and its taxonomic consequences in ferns. The two east¬ 
ern North American taxa, A. angustum and A. asplenioides. long have been 
perceived as close relatives separable by distinctive characters that are con¬ 
sistent within the vast northern and southern areas they respectively occupy. 
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Table 7 . Matrix of pairwise values for Rogers ( 1972 ) genetic similarity (above diagonal) and Nei 
( 1978 ) unbiased genetic identity (below diagonal). 


Population 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1. Mt. Ste. Hilaire, QE 

2. Barnet, VT 

3. North Hudson, NY 

4. Schenevus, NY 


5. 


0.975 0.962 0.953 0.939 0.853 0.865 0.847 0.862 0.881 


1.000 


0.956 0.938 0.930 0.851 0.864 0.853 0.858 0.879 


0.999 0.997 


0.944 0.943 0.844 0.856 0.840 0.852 0.873 


1.000 0.996 0.995 


0.942 0.841 0.842 0.831 0.843 0.863 


0.993 0.990 0.994 0.993 


0.817 0.821 0.803 0.825 0.843 


0.915 0.913 0.906 0.900 0.883 


0.924 0.921 0.916 0.908 0.892 0.996 


Ralston, PA 

6. Shirley, NJ 

7. Mountjoy Store, VA 

8. Hopewell, VA 

9 . Sandy Run Swamp, VA 0.924 0.921 0.916 0.905 0.891 0.998 1.000 0.998 
10. Pond Drain, VA 


0.943 0.951 0.959 0.932 


0.912 0.911 0.903 0.893 0.875 0.998 0.997 


0.950 0.966 0.959 

0.951 0.937 


0.956 


0.938 0.936 0.934 0.923 0.911 0.990 0.998 0.990 0.996 


but that intergrade and recombine to form a hybrid zone in their relatively 
narrow region ot overlap. This perception is supported by the data from 
the present study. The spores of the two taxa show striking and consis¬ 
tent differences in the perispore sculpturing, a low papillate perispore in 
A. cingustum versus a rugose perispore in A. asplenioides, and frequencies 
ot allozvmes exhibit strong difterences between as compared to within the 
taxa. An abrupt shift in allele frequencies at the four most polymorphic loci 
[Idh-t, Pgm-2, Tpi-2, and Pgi-2 ) ot Athyrium corresponds to the geographic 
boundary between A. cingustum and A. asplenioides, i.e., between northern 
Pennsylvania and southern New Jersey in our sample. Although virtually all 
alleles were shared between the two taxa, some alleles that predominated or 
were trequent in one taxon were nearly absent in the other, e.g., ldh-l A of 
A. asplenioides, Idh-l c ot A. angustum, and Pgi-2 Cl of A. angustum. In other 
cases, alleles were more extensively shared between the taxa but at very 
different frequencies, e.g. Idh-l 1 ', Pgm-2 B , Pgm-2 C , and Tpi-2 B (Table 3). 
UPGMA analysis ot Athyrium resulted in two distinct taxon clusters joined 
at a substantially lower similarity value (S = 0.849) than that joining popula¬ 
tions within their respective taxon clusters (S = 0.938 for A. angustum', 


Table 8 . Moans of pairwise values for Rogers’ Similarity (S) and Nei s Genetic Identity (I) tor 

comparison within and between Athyrium taxa. Ranges are given in parentheses. Each category 
was represented by ten pairwise comparisons. 


Taxon combination 

i 

S 

angustum-angustum 

0.996 

0.948 


(0.990-1.000) 

(0.930-0.975) 

asplen ioides-asp len i oi des 

0.996 

0.951 


(0.990-1.000) 

(0.932-0.966) 

angustum-asplenioides 

0.911 

0.848 


(0.875-0.938) 

(0.803-0.881) 
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Mt. Ste. Hilaire, QUE 
Barnet. VT 
North Hudson, NY 
Schenevus, NY 
Ralston-1, PA 
Shirley, NJ 
Hopewell., VA 
Mountjoy Store, VA 
Sandy Run Swamp NC 
Pond Drain, VA 



ang us turn 


asplenioides 


Rogers Similarity 


Fk;. 5. Dendrogram resulting 
laxity (Table 7). 


from UPGMA analysis based on pairwise values of Rogers’ Simi- 


S = 0.947, for A. asplenioides). These spore and allozyme differences for the 
most part are consistent among populations of each taxon, yet there is evi¬ 
dence of introgressive hybridization in their region of overlap as discussed 
below. 

The distinctness of these taxa and the consistencv of the character differ- 

•/ 

ences specified by Butters (1917) frequently have been questioned bv state¬ 
ments indicating that the taxa intergrade extensively (e.g., Benedict, 1934: 
Weatherbv, 1936: Fernald, 1946; Shaver, 1954; VVherrv, 1961). However, 
such statements seem anecdotal in that they are not accompanied by docu¬ 
mentation of character combinations in specimens. The most extensive 
specimen-based data set is that of Liew (1971, 1972), who carried out a phe- 
netic analysis of North American Athyrium sensu lato based on 170 speci¬ 
mens scored for 99 characters. Liew indicated that cluster analysis based on 

v 

paired similarity coefficients separated each Athyrium taxon, including 
A. angustum and A. asplenioides, into its own cluster, but the summary 
dendrograms presented leave it uncertain as to whether some of the speci¬ 
men placements were inappropriate or equivocal. 

In the present study, the degree of differentiation between A. angustum 
and A. asplenioides may appear exaggerated by the omission of localities 
where the two taxa co-occur. This omission was neither intentional nor an 
oversight, rather it resulted from a failure to discover such localities despite 
a substantial effort to do so. A search for co-occurring populations in south¬ 
ern Pennsylvania resulted in finding only a few isolated individuals of 
A. asplenioides in this highly agriculturalized and urbanized region, and 
no sizable populations from which allele frequency data could be obtained. 
All individuals in populations sampled for the present investigation from 
Quebec south to northern Pennsylvania were readily assignable to A. angus¬ 
tum on the basis of leaf and spore morphology. Similarly, all individuals in 
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populations from North Carolina north to northern Virginia were assignable 
to A. asplenioides. However, a genetic influence from A. angustum is sug¬ 
gested by the high frequency of characteristically A. angustum alleles Idh-1 8 

& ~ . r ™ _ _ 

and Pgi-2 ’ in the northernmost A. asplenioides population Shirley, NJ, as 
well as Pgm-2 B and Tpi-2 fi in the highest elevation A. asplenioides popula¬ 
tion Pond Drain, VA. Evidence of introgressive hybridization is augmented by 
the intermediate spore morphologies encountered at the Shirley, NJ, locality. 
No individuals assignable to A. angustum were found at this site, but it is 
possible that A. angustum spores coidd have migrated from populations fur¬ 
ther north and effected hybridization (Wagner, 1943). The prevalence of fully 
intermediate spores in six of the individuals, five of which are heterozygous 
for taxon marker allozymes, suggest that these are first-generation hybrids. 
The skewed array ot spore morphologies of other individuals suggests that 
they are backcrosses and provides evidence that hybridization has gone 
beyond the first generation resulting in introgression. Thus, spore and iso¬ 
zyme data combine to support previous morphology-based suggestions that 
the region ol overlap between the taxa, and possibly higher elevations in the 
southern Appalachians as well, could represent a hybrid zone resulting from 
secondary contact between recently diverged sister taxa (Benedict, 1934; 
Shaver, 1954; Sciarretta et al., ms). 

Narrow hybrid zones in which formation of fertile hybrids and backcrosses 
result in intergradation between divergent taxa occur between numerous spe¬ 
cies or subspecies pairs of animals and angiosperms (reviewed by Arnold, 
1997), and occur in a few agamosporous ferns (Gastony and Windham, 
1989). Hybridization between fern species usually results in spore abortion 
due to abnormal meiosis, and introgressive hybridization between homo- 
ploid taxa as divergent as A. angustum and A. asplenioides is decidedly rare 
in ferns, with only three cases having been documented previously: (1) 
swarms ot fully fertile hybrids involving Pteris cjuadriaurita and P. multia- 
urita, first generation as well as backcrosses, occur in disturbed forests of 
Sri-Lanka (Walker, 1958); (2) extensive hybridization among three species of 
the tree tern genus Alsophila resulted in a complex swarm ot fertile hybrids 
in 1 uerto Rico, a scenario hypothesized to give rise to new species through 
allogamous allohomoploidv (Conant and Cooper-Driver, 1980: Conant, 1990): 
(8) morphological and molecular data provide evidence ot hybrid swarms 
between PoWstichum munitum and P. inibriccms ot northwestern North 
America (Mayer and Mesler, 1993: Mulleniex et al ., 1998). Ot these, the situa¬ 
tion in Polystichum most closely parallels that ol the Athyriuni taxa anyus- 
tum and asplenioides . The two PoI\ r stichuni taxa share alleles at most 
enzyme loci, but are differentiated with respect to frequencies at three loci 
resulting in a lower magnitude of genetic identity between the taxa (1 = 

O. 842) than within each taxon (means: 1 = 0.974 for P. imhricans , 1 = 0.957 for 

P. munition; Soltis et a/., 1990, 1991). Additionally, the Polystichum taxa hy¬ 
bridize introgressively, although thev do not form a hybrid zone as in Athv- 

- ^ ^ 

rium, rather hybridization occurs at various localities across a broad area ot 
sympatry between the two taxa (Mayer and Mesler, 1993; Mulleniex et al.. 
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1998). Maintenance of distinction between these two Polvstichum species 
most likely results from a combination of their partial intersterility and di¬ 
versifying selection imposed by the very different habitats—forests versus 
exposed cliffs—occupied by P. munitum and P. imbricans respectivelv. 

The nature, frequency, and geographic extent of hybridization between 
A. angustum and A. cisplenioides remain uncertain and merit further 
research that combines field, herbarium, molecular, and breeding studies. It 
is critical to determine with greater precision and full documentation the 
degree to which these two taxa maintain versus blend their macromorpho- 
logical, micromorphological, and allozvmic differences where they coexist. It 
is unknown whether there is preferential mating within taxa. whether taxon 
characters tend to remain associated in the face of hybridization or are com¬ 
pletely recombined, and whether there exists a cline for allozvme frequen¬ 
cies and morphological characters within the overlap region. There is a need 
for intensive exploration for coexisting A. angustum and A. asplenioides pop¬ 
ulations in the areas of their overlap from the eastern seaboard through the 
midwest, and in the southern Appalachians where the existence of cryptic 
taxa has been hypothesized (Wagner and Wagner. 1966). Isozyme studies of 
these populations should be coordinated with critical analyses of morphologi¬ 
cal character combinations obtained from numerous specimens and with 
experimental crosses that can quantify the propensity of the taxa to hybridize. 

Beyond the taxonomic significance of hybridization, the unprecedented 
formation of normal intermediates between spores of such divergent mor¬ 
phology provides an opportunity to gain insight into the genetics underlying 
spore morphology (Schneller, 1989). The variability of spore morphology 
within individuals of the putative primary hybrids from the Shirley. N) site, 
which includes expression of parental types, indicates that inheritance is 
polygenic rather than a simple one-or-two-gene inheritance mechanism, and 
that the spore genotype determines or at least influences perispore pheno¬ 
type. Contrasting observations and inferences were obtained by Schneller 
(1989), who reported the formation of normal intermediate spores in experi¬ 
mentally produced hybrids between A. angustum and A. asplenioides . but 
found that all spores from a sporangium were of the same type, implicating 
sporophytic determination of the perispore morphology. Explanations that 
can account for these differing observations include the possibility that the 
Shirley, NJ, plants were not true FI hybrids, or that there may be a maternal 
effect that varies as a polymorphism. Clearly, further studies of the inheri¬ 
tance of spore morphology are in order. 

Taxonomic Conclusion: At What Rank Should A. angustum and A. asple¬ 
nioides Be Recognized ?—Over the second half of this centurv, the nature of 

w 

pteridophyte species has been clarified significantly by application of tech¬ 
nological advances in cytology and molecular systematics in combination 
with detailed morphological and field studies (Manton, 1950: Wagner, 1963: 
Haufier, 1987, 1989; Paris et al., 1989: Conant, 1990). Nonetheless, definitive 
ranking of closely related divergent taxa with overlapping ranges, such as 
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the two Athyrium taxa considered here, remains challenging due to the un- 


Mayr 


DeQueroz 


unpredictable nature of reproductive interactions between diverged taxa ex¬ 
periencing secondary contact (Arnold, 1997). The rank assigned to A. angus¬ 
tum and A. cisplenioides has varied considerably in floristic treatments 
published in this century. While Small (1938) and Wherry (1948, 1961) fol¬ 
lowed Butters (1917) in treating these taxa as separate species, other authors 
have tended to treat them as infraspecific taxa, either subspecies (Lellinger, 
1985) or varieties (Fernald, 1950; Mickel, 1979; Codv and Britton, 1989; 
Gleason and Cronquist, 1991; Kato. 1993). 

Spore and isozyme data combined indicate that populations of these two 
taxa are significantly divergent, exhibiting greater differences than ordinarily 
encountered within single species of ferns thus far studied. However, the 
fertility ot hybrids (Schneller, 1989) provides a potential for introgression 
between the two taxa. The preliminary evidence that hybridization and 
introgression do in fact occur indicates that A. angustum and A. asple- 
nioides would be considered conspecific under species definitions as differ¬ 
ent as the Biological Species Concept (Mavr, 1942) and the Phylogenetic/ 
Diagnostic Species Concept (Cracraft, 1983; Davis and Nixon, 1992). None¬ 
theless. in practice numerous pairs of taxa that form hybrid swarms or zones 
are treated as distinct species (Grant, 1981; Arnold. 1997). 

The occurrence ot northern and southern infraspecific taxa in eastern 
North American Athyrium is paralleled in Pteridium aquilinum L., which 
comprises northern and southern varieties latiusculum (Desv.) Underw. and 
pseudocaudatum (Clute) Heller, respectively, and which shows north to 
south dines in allele frequencies near the overlapping taxon boundary 
(Speer et ol., 1998). However, the pattern ot allozyme variation in Pteridium 
(litters trom that in Athyrium in that the most abrupt shift in allele frequen¬ 
cies occurs within the range ot P. latiusculum, genetic identities between 
populations ot the two Pteridium varieties ranged substantially higher (1 = 
0.916 to 0.999) than those of the two Athyrium taxa (1 = 0.875—0.938. mean 
= 0.91 1), and UPGMA clustering tailed to separate the varieties (Speer et ah, 
1998). On the basis of the lack of genetic differentiation between P. latiuscu¬ 
lum and P. pseudocaudatum, variety was indicated as the highest rank at 
which to recognize them (Speer et ah, 1998). In contrast, the more substan¬ 
tial ditterentiation between the Athyrium taxa angustum and asplenioides 
and the consistent characters uniting a vast number ot individuals north and 
south ol their hybrid zone suggest that the taxa should be ranked at least at 
the level ol subspecies. Ranking at the species level would not be inconsis¬ 
tent with the treatment of such taxa in the broader plant literature. 
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